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ABSTRACT 

We show that S-minute acoustic oscillations may res- 
onantly convert into Alfven waves in the /3^1 region 
of the solar atmosphere. Considering the 5-minute 
oscillations as pumping standing acoustic waves os- 
cillating along unperturbed vertical magnetic field, 
we find on solving the ideal MHD equations that am- 
plitudes of Alfven waves with twice the period and 
wavelength of acoustic waves exponentially grow in 
time when the sound and Alfven speeds are equal, 
i.e. Cs ~ va. The region of the solar atmosphere 
where this equality takes place we call a swing layer. 
The amplified Alfven waves may easily pass through 
the chromosphere and transition region carrying the 
energy of p- modes into the corona. 

Key words: solar atmosphere; 5-minute oscillations; 
Alfven waves. 



1. INTRODUCTION 



It is generally considered that solar 5-minute pho- 
tospheric acoustic oscillations cannot penetrate into 
the upper regions due to the acoustic cutoff of strati- 
fied atmosphere. For the typical photospheric sound 
speed Cs = 7.5 km/s the cutoff frequenc y is 0.03 
s~^, which gives the cutoff period of 210 s l|R,obert j . 
I200I . This means that the sound waves with 5 
min (~ 300 s) period are evanescent. Acoustic os- 
cillations cannot penetrate into the corona also due 
to the sharp temperature gradient in the transition 
region. On the other hand, 3 and 5-minutc in- 
tensity oscillations are intensively observed in the 
corona by the space satellites SOHO (Solar and He- 
liospheric Observatory) and TRACE (Transitio n Re- 
gion and Coronal Explorer) l|De Moorte]ll2002^ . Re- 



cently 'De Pontieii et all (J2005') have discussed how 
photospheric oscillations can be channelled into the 
corona through inclined magnetic fields. Another so- 
lution of this controversy is the conversion of acoustic 
oscillations into another wave mode, which may pass 
through chromosphere/transition region. 

Recent two dimensio n al numerical simul ations 
(|R,osenthal et alJ . l2002t iBogdan et all |2Q03lf) out- 
lined the importance of /3~1 region in the solar at- 
mosphere. They found the coupling of fast and slow 
magnetosonic waves at this particular region. Re- 
cent modelling of the plasma f3 in the solar atmo- 
sphere (Gary 2001, see Fig. 3 of that paper) shows 
that VA ^ Cs, i.e. /3 ~ 1 (actually /3 ~ 1.2 for 
7 = 5/3), may takes place not only in lower chro- 
mosphere, but also at relatively low coronal heights 
(e.g., ~ 1.2i?o from the surface, whe re Rn is the so- 
lar ra dius). Also latest observations (Mugl ach et al.l 
l2005il suggest the possible transformation of com- 
pressible wave energy into incompressible waves at 
/3 « 1 region of the solar atmosphere. Thus this 
particular region may be of importance due to con- 
version of compressible wave energy into incompress- 
ible Alfven waves (or into MHD kink waves in thin 
photospheric magnetic tubes). 

The coupling of propagating sound and Alfven waves 
at /3~1 has been proposed by Zaqarashvili and 
Roberts (2002,2005). They found that a sound wave 
is nonlinearly coupled to the Alfven wave with dou- 
ble the period and wavelength when the sound and 
Alfven speeds are equal, i.e. Cs ~ va- 

lUlrichI l)1996() has reported observations of Alfven 
waves in the solar photosphere and lower chromo- 
sphere with substantial power at frequencies lower 
than the 5 minute oscillation. In the power spec- 
trum of magnetic oscillations (Fig. 3 in that paper) 
there is significant power at about 10 min. 



Here we show that standing acoustic waves oscil- 
lating along uniform magnetic field lines effectively 

generate Alfvcn waves with double their period and 
wavelength in the /3^1 regions of the solar atmo- 
sphere. The case of propagating waves is discussed 
in Zaqarashvili and Roberts (2002,2005). 



2. STATEMENT OF THE PROBLEM 
AND DEVELOPMENTS 



Consider fluid motions u in a magnetised medium 
(with zero viscosity and infinite electrical conductiv- 
ity), as described by the ideal MHD equations: 



— + (u-V)B = (B-V)u - BV-u, 



p— -h p(u-V)u ^ 



dt 
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where p is the medium density, p is the pressure, u is 
the velocity, B is the magnetic field and 7 is the ratio 
of specific heats. For simplicity the stratification is 
neglected here, but we plan to take it into account 
in future. 

Cartesian coordinate system is adopted with the z 
axis directed vertically upwards from the solar sur- 
face. Spatially inhomogeneous (along the x axis) 
magnetic field is directed along the z axis (see Fig.l), 
i.e. 

So = (0,0,B,(a;)). (5) 

Plasma pressure and density also are assumed to 
have X dependence, so they are expressed as po(a^) 
and pq{x) respectively. The magnetic field and pres- 
sure satisfy the transverse pressure balance condition 



Pq{x) + 



Stt 



const. 



(6) 
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V- Standing acoustic oscillations 
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Figure 1. Schematic picture of the proposed process. 
The z coordinate is directed vertically upwards, while 
the X and y axes are parallel to the solar surface. The 
m,agnetic field and density are inhomogeneous along 
the X coordinate. Acoustic waves propagate up and 
down, i.e. along the z axis. Alfven wave velocity 
oscillates along the y axis. Standing acoustic waves 
are resonantly converted into Alfven waves near the 
region where va ^ Cg called as swing layer. 



yz plane. Then only sound and Alfven waves arise. 
The velocity component of sound wave is polarised 

along the z axis and the velocity component of the 
Alfven wave is polarised along the y-axis. In this 
case equations (l)-(3) take form 
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Plasma /3 is defined as 



87rpo(a;) 



= 2c?/74(a^), 



(7) 



where Cg = \J jpo / po and VAix) = Bz/^4:TTpo are the 
sound and Alfven speeds respectively. Note that we 
suggest the temperature to be homogeneous so the 
sound speed does no depend on the x coordinate. 

We consider wave propagation along the z axes (thus 
along the magnetic field) and wave polarisation in 



where p = Po + Pi and p = po + Pi denote the to- 
tal (unperturbed plus perturbed) pressure and den- 
sity, Uy and Uz are the velocity perturbations (of the 
Alfven and sound waves, respectively), and by is the 
perturbation in the magnetic field. Note, that in 
these equations the x coordinate stands as a param- 
eter. 

Acoustic waves oscillate along the z axis, so that the 
velocity component has nodes at points z = and 
z = I, i.e. Uz = at z = 0,1, and thus we take 



Uz = v{t) sm{ksz), 



(13) 
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pi = pi{t)cos{ksz), (14) 
where kg is the wavenumber of sound wave such that 

ksl = = riTT, 
As 

so 

•As 

where n — 1, 2.... 

We express the Alfven wave components as 

by = b{t) cos{kAz), (15) 

Uy = u{t) siii{kAz), (16) 
where kA is the wavenumber of the Alfven waves. 

Substitution of expressions (13)-(16) into equations 
(8)-(12) and averaging with z over the distance (0, 1) 
leads to the cancelling of all nonlinear terms, which 
means that waves do not interact. However in par- 
ticular case, when wave numbers kg and kA satisfy 
the conditions 

ks = '2k A, 



equations (8)-(12) take the form: 

8h k A 

1^ = kAB^u{t) - -fv{t)b{t). 
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which means that the waves may interact as the non- 
linear terms remain. 



Equation (23) reflects the fact that the Alfven speed 
is modified due to the density variation of stand- 
ing acoustic wa ve. The similar equ ation for (3^1 
was obtained bv lZaaarashvilil l)2001|) . It is seen from 
this equation that the particular time dependence of 
density perturbation determines the type of equation 
and consequently its solutions. If we consider the 
initial amplitude of Alfven waves smaller than the 
amplitude of acoustic waves, then the term with 
in equation (20) can be neglected. This means that 
the backreaction of Alfven waves due to the pondero- 
motive force is small. Then the solution of equations 
(20)-(21) is just harmonic function of time 



pi = apocos{ujst), 



(24) 



where uJs is the frequency of standing acoustic wave 
and a <C 1 is the relative amplitude. Here we con- 
sider the small amplitude acoustic waves a <C 1, so 
the nonlinear steepening due to the generation of 
higher harmonics is negligible. Then the substitu- 
tion of expression (24) into equation (23) leads to 
the Mathieu equation 



d^b ,2 2 
dt^ + ^^^^ 



a ■ 
— cos(wsi) 



0. 



(25) 



The solution of this equation with frequency 
ujs/2 has an exponentially growing charac- 
ter, thus the main resonant solution oc- 

I200I 



curs when l|Zaq arashvil i and Robert^ 
IShergelashvih et all l2005il 



UA 



VAkA = 



(26) 



where uja is the frequency of Alfven waves. Since 
kA = ks/2, resonance takes place when 



VA 



(27) 



Substitution of u from equation (18) into equation 
(19) and neglecting of all third order terms leads to 
the second order differential equation 



kAV 



db 
di 



k\v\ 



2/30 



Pi 



kA dv 
~2"dt 



= 0. (22) 

This equation describes the time evolution of Alfven 
wave spatial Fourier harmonics expressed by (15)- 
(16) forced by standing acoustic waves. 

In this equation, the first derivative with time can 
be avoided by substitution of function 

b = b{t)e-^I'"'\ 
which after dropping third order terms leads to 



k\v\ 



2po 



Pi 



0. 



(23) 



Since the Alfven speed va{x) = Bz{x)/ y/ 47rpo (x) is 
a function of the x coordinate, then this relation is 
satisfied at a particular location along the x axis (see 
Fig. 1). Therefore near this region the acoustic os- 
cillations will be resonantly transformed into Alfven 
waves. We call this region the swing layer, by anal- 
ogy with mecha nical swing interactions (see a similar 
consideration in lShergelashvili et alJ l)2005() '). 

Under condition (26) the solution of equation (25) is 



b{t) = 6oe^* 



UJs . UJs ' 

COS — t=Fsm — t 
2 2 . 



(28) 



where bo = b{0) and the phase sign depends on a; it 
is -|- for negative a and — for positive a. 

Note that the solution has a resonant character 
within the frequency interval 



(29) 
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Figure 2. Numerical simulations of wave conversion. 
Upper panels show the density and velocity compo- 
nents of standing acoustic oscillations, while lower 
panels show the magnetic field and velocity compo- 
nents of Alfven waves. Alfven waves with twice the 
period of acoustic oscillations grow in time. 

This expression can be rewritten as 



Thus the thickness of the swing layer depends on the 

acoustic wave amplitude! 

Therefore the acoustic oscillations are converted into 
Alfven waves not only at the surface va = Cg but also 
near that region, namely at 

VA = c, (l ± ^) . (31) 

Thus the resonant layer can be significantly wider for 
stronger amplitude acoustic oscillations. 

Numerical solution of equations (18)-(21) (here the 
backreaction of Alfven waves is again neglected) is 
presented on Fig. 2. The amplification of Alfven 
waves with double the period of acoustic oscillations 
is clearly seen. 



4. ACKNOWLEDGEMENTS 

D. Kuridze and T. Zaqarashvili acknowledge the fi- 
nancial support from conference organisers. 



REFERENCES 

Bogdan, T.J., Hansteen, M., Carlsson, V. et al. 2003, 

ApJ 599, 626 

De Moortel, I. Ireland, J., Hood, A. W. and Walsh, 
R.W. 2002, A&A 387, L13 

Do Ponticu, B., Erdelyi, R. and De Moortel, I. 2005, 

Ap.J, 624, L61 

Gary G.A. 2001, Solar Phys. 203, 71 

Muglach, K., Hofmann, A. and Staude, J. 2005, 

A&A, 437, 1055 

Roberts, B. 2004, In Proc. SOHO 13 'Waves, Os- 
cillations and Small-Scale Transient Events in the 
Solar Atmosphere: A Joint View from SOHO and 
TRACE', Palma de Mallorca, Spain, (ESA SP- 
547), 1 

Rosenthal, C. S., Bogdan, T. J., Carlsson, M., Dorch, 
S. B. F., Hansteen, V., Mcintosh, S. w'., McMurry, 
A., Nordlund, . and Stein, R. F. 2002, ApJ 564, 
508 

Shergelashvili, B.M., Zaqarashvili, T.V., Poedts, S. 
and Roberts, B. 2005, A&A 429, 767 

Ulrich, R.K. 1996, ApJ. 465, 436 

Zaqarashvih, T.V., 2001, ApJ, 552, LSI 

Zaqarashvih, T.V. and Roberts, B. 2002, Proc. 10th. 
European Solar Physics Meeting, "Solar Variabil- 
ity: ^From Core to Outer Frontiers" . Prague, 
Czech Republic, 9-14 September 2002 (ESA SP- 
506, December 2002), 79 

Zaqarashvili, T.V. and Roberts, B. 2005, submitted 



3. CONCLUSIONS 



We suggest that 3 and 5-minute acoustic oscillations 

in photosphere/chromosphere can be resonantly con- 
verted into Alfven waves, or possibly into MHD kink 
waves in thin photospheric magnetic tubes, this pro- 
cess acting in the region of the solar atmosphere 
where VA^Vg. Generated transversal waves may 
then propagate through the transition region into 
the corona, where they can deposit their energy back 
into density perturbations. The process can thus be 
of importance in coronal heating. 



